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Green Operators of Networks with a new vertex 

A. Carmona, A.M. Encinas, S. Gago, M. Mitjana 

Abstract. Any elliptic operator defines an automorphism on the orthogonal 
subspace to the eigenfunctions associated with the lowest eigenvalue, whose 
inverse is the orthogonal Green operator. In this study, we show that elliptic 
Schrodinger operators on networks that have been obtained by adding a new 
vertex to a given network, can be seen as perturbations of the Schrodinger 
operators on the initial network. Therefore, the Green function on the new 
network can be computed in terms of the Green function of the original 
network. 


1. Introduction 

Discrete elliptic operators can be seen as the discrete counter part of elliptic partial differential opera¬ 
tors. In particular, positive semi-definite Schrodinger operators defined on a finite network are examples of 
those self-adjoint operators. Any elliptic operator defines an automorphism on the orthogonal subspace to 
the eigenfunctions associated with the lowest eigenvalue, whose inverse is the orthogonal Green operator. 

In [4], some of the authors analyzed the effect of a perturbation of the network by computing the 
effective resistance of the perturbed networks through Sherman-Morrison-Woodbury like-formulas, instead 
of using the Sherman-Morrison formula recursively. In fact, since adding edges to a network does not 
modify the space of functions on the vertex set of the network, this class of perturbation was placed into the 
general framework of perturbations of discrete elliptic operators. Specifically, we showed that this problem 
corresponds with the superposition of rank one perturbations that are orthogonal to the eigenfunction 
associated with the lowest eigenvalue of the elliptic operator. 

The scenario changes when the perturbation consists on adding new vertices to the network. Only 
few works have tackled the problem of adding a new vertex, see for instance[2]. In this work, we consider 
perturbations that consist on adding a new vertex to a network. After some well-known operations on the 
Schrodinger operator of the perturbed network, that involves the inverse of the Schur complement of the block 
corresponding to the added vertices, we show that this Schur complement can be seen as a perturbation 
of the Schrodinger operator of the original network, understood as a discrete elliptic operator, that is a 
superposition of rank one perturbations that, this time, are not orthogonal to the eigenfunction associated 
with the lowest eigenvalue of the elliptic operator. Therefore, we can apply the general theory developed in 
[4] for this kind of perturbations. 

We start the study by revisiting the perturbation of an elliptic operator with a sum of projections that 
can be, or not, orthogonal to the eigenfunction associated with the smallest eigenvalue. Thus, we consider 
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the relation between the Green operator of the new operator in terms of the Green operator of the previous 
one. 

Next section is devoted to the application of the mentioned results to the addition of a new vertex to the 
network F in order to get a network Fb Moreover, we obtain the relation between the Schrodinger operators 
of the two networks F and F' and in addition, we give the explicit expression of the matrix associated with 
the Green operator. 


2. Specific notation and preliminary results 


Given a finite set V of n elements, we denote by C(V) the space of real valued functions on V. For any 
function u G C(V), the associated vector in K" will be denoted by u. For any vertex x € V, the Dirac function 
at X is denoted by € C(V); the scalar product on C(V) is (u,v} = for each u,v G C(V). A 

unitary and positive function w is called a weight and D(D) denote the set of weights. 

If JC is an endomorphism of C(V), it is self-adjoint when {IC{u),v) = {u,IC{v)) for any u G C{V). 
Moreover, /C is positive semi-definite when (/C(m),u) > 0 for any u G C{V). A self-adjoint operator /C is 
elliptic if it is positive semi-definite and its lowest eigenvalue A is simple. Moreover, there exists a unique 
unitary function uj G C(V), up to sign, satisfying IC{oj) = Xoj, so JC is called (A,a;)-elliptic operator. It is 
straighforward that a (A, a;)-elliptic operator is singular iff A = 0. We denote by A'^ = A“^ iff A ^ 0, and 
A1 = 0 otherwise. 

Any function AT : D x G —?> K is called a kernel on V, and it determines an endomorphim of C(V) 
by assigning to any u G C{V) the function lC(u) = Eyeu Conversely, each endomorphism of 

C{V) is determined by the kernel K(x,y) = {IC{ey),ex) for any x, y gV. Therefore, an endomorphism JC is 
self-adjoint iff its kernel K is a symmetric function. 

Given a, t G C(V), we denote by Vcr,T the endomorphism of C(V) that assigns to each u G C{V) the 
function Va^riu) = {t,u)u, and it is called a projector, as it assigns to any u G C(V) its projection on a 
along T. Observe that the corresponding kernel is Pa,T{x,y) = {u ® T){x,y) = a{x)T[y). In particular, when 
a; ^ 0 the endomorfism 'Pu),u: is denoted simply by P^i- 

Given A > 0, a; € fl(fo) and a (A, a;)-elliptic operator P, we will be concerned with the so-called 
Poisson equation for P on V: for a given / G C{V) find u G C(V) such that P(u) = f. As P defines an 
automorphism on oj^, the inverse of a (A, a;)-elliptic operator P on w*- is called orthogonal Green operator 
and it is denoted by G. This operator on uj^ can be extended to CiV) by assigning to any / G C(V) the 
unique solution of the Poisson equation P{u) = / — Vujif)- 

Now consider a non-null function a G C(V), the associated self-adjoint projection P^r and the operator 
Ho- = 'Per, called perturbation of P by cr. The relation between Green operators of both Per and P can 
be found in [H Corollary 3.6]. 


Corollary 1. Consider a G C{V). If X = (cr, uj) = 0, tJien 

1 


G'^’’ = G - 




1 + (S(cr),cr 

wJiereas wJien eitJier X > 0 or {a, uj) ^ 0, tJien is invertible and 

1 r 


n 


-1 


— ^ ~ ~ {PgP),^ ~'Pur,g{a)) - {l+ {Q{(^)p))'J-’u. 


where /? = A(1 -1- (^(cr),cr)) -|- (cr, w)^. 


Moreover, if we consider ai G C(V), i = 1,..., m + £ such that tJi ^ uj'‘~ for i = 1,... ,m and rim+i C w*", 
i = 1,... ,£ the operator 

n = p+'^Pa., 


i=l 

is a perturbed operator. Here, m or i can be equal to 0. The relation between the corresponding inverse 
operators is given in the following theorem. 
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Theorem 1. [H Theorem 3.5] The operator TL is positive semi-definite, positive definite when m > 1 and 
moreover 

m-\-£ m-\-£ 

= G + hVuj + ^ ^ f^ijrg(cri),g{crj)j 

i—1 

where [bij) = (/ + 0 -^))“^, and 


\ r,s—l / 

m 

hi = h bir{(Tr,Lu), i = 1,... ,m + £, 

r—1 

^ij ^ij ^ f ^ ^ ^ir j f ^ ^ ^jr {^r j^)j; j ^ T ^ 


\r—l 


\r—l 


On the other hand, the Schur complement provides us with a fundamental tool for the results on next 
sections. 

Lemma 1. If A £ A^nxn, B G TWnxm, D G TWmxm invertible, and S = A — then 

(A B y _ / St -StBD-i \ 

BT D j “ -D-^B^St D-i + D-iBTStBD^O ’ 

where stands for the Moore-Penrose inverse of the matrix C. 


3. Adding a new vertex 

In the following, the triple F = {V, E, c) denotes a finite network, i.e., a connected graph without loops 
nor multiple edges with vertex set V, with cardinality n, and edge set E, in which each edge Cxy G E has 
assigned a value c(x, y) > 0 named conductance. The conductance c is a symmetric function c : V x V —>■ 
[0, oo) such that c{x,x) = 0 for any x £ V and where the vertex x is adjacent to vertex y iff c(x,y) > 0. 
Given a weight on V, w G fI(V), for any pair of vertices (a:, y) G V x V the w-dipole between x and y is the 
function Txy = The Laplacian of the network F is the endomorphism of C{V) that assigns to 

each u £C{V) the function 

C{u){x) = '^c{x,y)[u{x)-u{y)], x £ V. 
y&v 

The Laplacian is a singular elliptic operator on C{V) and moreover C{u) = 0 iff u is a constant function. 
Given q £ C(V) the Schrddinger operator on F with potential q is the endomorphism of C{V) that assigns 
to each u £ C(V) the function £q{u) = £{u) + qu. Given a weight w G fl(V), the potential determined 
by uj is the function q^^j = —^£{u}). It is well-known that the Schrodinger operator £q is (A, w)-elliptic iff 
Q = Qui + A, see [T]. Moreover, it is singular iff A = 0 and then, £q^{v) = 0 iff u = aw, o G K. We denote by 
G\,u] the orthogonal Green operator associated with £q and by G;^ its corresponding kernel. From now on, 
we consider fixed the value A > 0, the weight w G r2(V) and the Schrodinger operator £q with q = q^^ + A. 

In this study we worry about perturbations of £q performed by adding a new vertex. Namely, let x' be 
a new vertex and assume we connect it to m vertices Xi,..., Xm G 1^, where 1 < m < n. The new network 
F' = {V, E', c') has vertex set V = V U {x'}, edge set E' = EU {ox^x', • ■ • > and conductance 

{ Oi > 0 if X = Xi £ V, y = x', 1 < i < m, 

c(x, y) if x,y £ V, 

0 otherwise. 

If uj(x') is a positive value assigned to x', we define a weight on F' by uj'(x) = U!(x)/^1 -j- U!(x')^, for 
any x £ V'. Observe that in this case it holds w(a;)/w(y) = w'(a;)/w'(y), for any x, y £ V'. 
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The following notation is useful in what follows. For any i = 1,..., m, we denote by pi = aiUj{xi)oj{x'), 


O'* = 


Pi \-^m ^ , f \-^m 9 

cr = and a = A + Y.i=i Pi- 


j{Xi) 


r'K) 


Proposition 1. If CJ is the Laplacian ofT' and p = Qu,' + A, where q^ji = -, then 

i. / 


UJ' 


C'p = Cq+Y. 'P^r, - Va,s^, On V, 


P'p ~ ~Ps^i,a + 




Proof. For any u G C(V') it is satisfied that 

C'u(xi) = Lu(xi) + ai{u{xi) — u(x')), i = 1,..., m, 


L'u{x') = Y ai(u{x') - u{xi)), 

i—1 

L'u{x) = Cu(x), otherwise, 


and in particular, 


qu,'{xf) = qu,{xi) - Oi + ( 


K^i) 


, ■j=l,...,TO, 


quj'ix') = -Yai-\ — rn E ai<-^ixi), 
z=l U}{x'} 

qu)'{x) = quj{x), otherwise. 


Therefore, 


Pi 


C u{xi) = Cqu{xi) H- Y~~^'^i^i) ~ 0 ,iU{x'), i = 1,..., m, 

^ UJ"^[Xi) 

, m m 

£pu{x') = Xu{x') + 3 

£pu(x) = Cqu{x), otherwise. 


and the result follows. 


The relation between the matrices associated with both Schrodinger operators of F and F' is given by 


l; = 


H -s 


-<;T 


S' a 


where H is the matrix associated with the operator 

m 

n = Cq + Y,P<^.^ 

and s is the column vector s = YiLi ('-i^xi, where for z = 1, ..., n are the vectors of the canonical basis. 

In order to compute the Moore-Penrose inverse of we will use Lemma and it will be useful to 
introduce the following perturbations. We define, ior k = 1,... ,m, 

A 


'^k — \ ^k') 

a 


^2/77,_1_zh 

and for i = 1,..., TO — 1, j = i + 1,..., m, let /c =---h j and 


'^k — 


1 PiPj ( ^Xi 


y/auj(x') \Uj(Xi) Uj{xj) J 
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(l;)^ = 


Theorem 2. The Moore-Penrose inverse of is given by 

M -iMs \ 

-isTM i + 4sTMs / ’ 

where M is the matrix associated with the operator 

m(m+l) rTx(m+l) 

2 2 

G P hP^ ^ ^ (-TTi),0; '^a;,0(7ri)] ^ ^ 

i=l i,j=l 

and where if {bij) = (J + (C/(7rj), 7ri))“^, 


h = X^a a + ^ brsPrPs 

\ r.s=l / 

A “ 


hi — ^ ^ ^irPr^ i — 1? ■ ■ ■ ; 


m{m + 1) 


h — h i h ] f h 1 ■ ■ — 1 Tn{rn + 1 ) 

hij — bij I / ^ birPr j I / , bjsPs I ) IjJ — 1) • • ■ ) „ 


\r—l 


\s^l 


Proof. From Lemma we have that 

(LP^ = 


St 


-Sts 

a 


-sTSt A + :^sTSts 


where S = H-s(8)sis the matrix associated with the operator 

a 


S = Cg + yP,^ - -Va- 
_ a 


Now, let us prove that 




f27Ti _ 1 _ 1 j \'1 

where aij = ^/a^rk^ k = ^ + j for i = 1,..., to — 1, j = z + 1,..., to. If Po- denotes the kernel 

of Va and P denotes the kernel of the operator in the right side of the equality, the claim is equivalent to 
prove that P — P^. 

Since P^ = a ® a, we have 

m 

Per 'y ^ ^ ^Xj)- 

hi=i 

On the other hand, for any z = 1,..., to, we have 

Per. = , ... (Sx, 0SxJ = - j -P (e^i 0£xj, 


oj(xiy 


Uj(Xi) 


and hence, 


V p., = u;(x') V -y (£x, 0£xj. 
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Moreover, for 1 < i < j < m, 

Pa.. = 


2 2 2 2 

{£xi ‘^£xi)+ , ..o , —^ {£xj ® £x 


Uj{x'Yuj{XiY ' ' Uj(x'Yuj(XjY 

{£xi <S>£xj +£xj <S>£xi) 


2 2 

pfPj 


Uj{x'Yuj{Xi)uj{Xj) 

a^ajUj{xj) aiajUj{xi) 

, \ \£xi Qv £xi> + /X 

Uj{Xi) 

0.iClj(^£xi ® £xj P £xj ^£xi^f 


and hence, 




aiajUj{xj) 


, , r V~'V~' OiiUjivlXi) , , 

{£xi ® Ea;;) + 2 ^ ^ w(ar) ^ 

5j 


l<2<j<m 


l<i<j'<m 


— ^ ^ aja;(Xj)^ (^a;. 0 ^a^i) CliCljiSxi ® Sxj) • 


3^i 


^J = l 


Taking into account that 

m m 

ajUj{xj) = ajUj(xj) — aiUj{xi) = {a — A)uj(x') — aiUj(xi), 


j=i 


i=i 


we obtain that 


Therefore, 


y ' Po-ij — (o; A) y ^ Pai y ^ aiaj (f 

l< 2 <j'<m 2=1 —1 


Xi 'ey ^x 


)■ 




5 = /:„ 


E 


fc=l 


Finally, from Theorem we get that = M and the result follows. 
Next we describe the coefficients of matrix {{G{T^e),T^k))- 

Lemma 2. The elements of the matrix A = {{G{Tre),TTk)) are given by 

! A\ ApkPl GiX u>{Xkj Xf) 

{A)k/ = : —7irT“7TT’ K,t=l,...,TO, 


a u{xk)uj{x() 


{PYk,i = 


^/APkPiPj 

aui{x') 


G>X^aj{Xkf xY Gix^ijj{Xki XY 
Lo{xk)oj{Xi) uj(xk)uj(xj) 


, „ (2m — 1 — i)i 

k = l,...,m, 1= ---h j, 


(a4)fcy = 


^X^uj(Xi^Xr( L^qj(x2,Xs) (jx^uj(Xj ^ xY ^ Gjx^u](Xj ^ X 


Pi Pj Pr Ps 

aUl(x'Y llx!{Xi)uj(Xr) Uj(Xi)uj(Xs) Uj(Xj)uj(Xr) Uj(Xj)uj(Xs) 

(2m — 1 — r)r 

forr = l,...,m — l, s = r + 1,... ,m, k = -^-h s, 

(2m — 1 — i)i 

and for i = I,... ,m — 1, j = i + 1,... ,m, 1= ---h j. 


□ 


2 
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Observe that we can deduce two special cases of Theorem if m = 1 that means the addition of a 
pendant vertex to the network (see 0), and if m = n it represents the join of a the new vertex with the 
graph; (see HD- 


Corollary 2. If m = 1 and a = c{x, x'), then 

(l;)' = ( 


M 


-xMe, 


- + 


where M is the matrix associated with the operator 

1 




AT’g(^) - Px {Vg(^a),u, - 'Pcj.Qicr)) - (1 + (a - A)GA,(^(a:, x)) Va. 
where ft. = A[1 + (a — A)G;^^^(x,a;)] + p^. 
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